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ABSTRACT 

We investigate the correlation between nine different dark matter halo properties using 
a rank correlation analysis and a Principal Component Analysis for a sample of haloes 
spanning five orders of magnitude in mass. We consider mass and dimensionless mea- 
sures of concentration, age, relaxedness, sphericity, triaxiality, substructure, spin, and 
environment, where the latter is defined in a way that makes it insensitive to mass. We 
find that concentration is the most fundamental property. Except for environment, all 
parameters are strongly correlated with concentration. Concentration, age, substruc- 
ture, mass, sphericity and relaxedness can be considered a single family of parameters, 
albeit with substantial scatter. In contrast, spin, environment, and triaxiality are more 
independent, although spin does correlate strongly with substructure and both spin 
and triaxiality correlate substantially with concentration. Although mass sets the scale 
of a halo, all other properties are more sensitive to concentration. 

Key words: cosmology: theory - cosmology: dark matter - galaxies: haloes 



> 

in 

cn 
o 



> 

X 
C3 



1 INTRODUCTION 

Hierarchical structure formation in the ACDM universe pre- 
dicts that galaxies form within virialised dark matter (DM) 
haloes which merge to form larger ones. Understanding the 
structure of dark matter haloes is therefore of crucial im- 
portance for models of the formation of galaxies. A large 
number of simulations have been used to study the correla- 
tions between halo properties like mass, concentration, ac- 
cretion history, spin, shape, substructure and environment 
(e.g. Navarro, Frenk, fc White] 



1997; 



Harkcr ct al. 2006; Bctt ct al. 2007 



Bullock et alj 2001 



Maccio et al.|2007 



IHahn 



et al ||2007| |Duffy et aJ.|2008||Wang et al.|2010| |Haas et al 



20111 and to investigate the dependence of halo clustering 



on several of these parameters (e.g. Gao, Springel, & White 


2005 


Wechsler et al.||2006| |Bett et al.||2007| |Gao & White 


2007 


Jing, Suto, & Mo|2007| | Wetzel et al.|2007||Li, Mo, & 


Gao 


2008 1 . Semi-analytic and halo occupation models usu- 



ally assume halo mass to be the most important or even the 
only property that defines the halo and its galaxy content 
(e.g.|Mo fc White||1996l |Cooray fc Sheth||2002[). However, 



there is mounting evidence that formation time is likely to be 
another very important parameter ( |Croton, Gao, fc White| 
2007||Taylor|2011[ ). 

In this work we present a numerical study of the rela- 
tions between halo mass, formation history, shape, dynam- 
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ical state, and environment at redshift z = 0. We make use 
of both a correlation and a Principal Component Analysis 
(PCA) for a set of 9 parameters. By combining 5 large, cos- 
mological N-body simulations, we are able to explore a mass 
range that spans 5 orders of magnitude with no particular 
mass dominating the distribution. 

This paper is organised as follows. In Section [2] we de- 
fine the halo properties that we consider and we discuss the 
simulations in Section[3] In Section|4]we discuss the analysis 
methods and present the results of our investigation. 



2 DEFINITIONS OF HALO PROPERTIES 

We will consider the 9 halo properties described below. Ex- 
cept for mass, all are dimensionless. 

Virial mass, M200 The virial mass is the mass within a 
sphere of radius -R200 centered on the most bound particle 
of the halo, where -R200 is the radius within which the over- 
density is 200 times the critical density of the Universe. 

Concentration, C200 The concentration of a halo is 
found by fitting an NFW density profile ( jNavarro, Frenk,] 
fc White]|1997 l to the particle distribution and calculating 
C200 = -R200 /r s where r s is the NFW scale radius. The values 
were obtained from puffy et al.|p008| . 

Age, zo.5 The age of a halo is defined as the redshift at 
which at least half of the final mass, i.e. 0.5 M200, is in FOF 
haloes with at least one-tenth of the final mass. The redshift 
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and mass fraction are calculated by linear interpolation be- 
tween adjacent snapshots. 

Relaxedness, r The relaxedness parameter is defined as 
the distance between the most bound particle and the centre 
of mass of the FOF halo divided by i?200- A more relaxed 
halo has a smaller r. 

Sphericity, s The halo sphericity, s = c/a, is denned as 
the ratio of the length of the minor axis c to the length of 
the major axis a. The axis lengths are calculated by diago- 



(2007) 



nalising the inertia tensor defined as in Bett et al. 

Triaxiality, T Triaxiality, T = (a 2 - b^)/(a z - c 2 ), where 
b is the length of the intermediate axis of the diagonalised 
inertia tensor. T — > 1 for a highly prolate halo and T — > 
for a highly oblate halo. 

Substructure, / sub The substructure parameter, / su b = 
M su b/M2oo, is defined as the fraction of the mass within 
^200 that is in substructures. We consider only sub-haloes 
with maximum circular velocities larger than 10 -0 ' 5 times 
that of the parent halo. 

Spin, A We use the modified spin parameter, A = 
j(\/2V2oo^20o) _1 , where j is the specific angular momen- 
tum and V200 the circular velocity of the halo ( Bullock et 
|al.|200l| . For details and a comparison of different spin pa- 



rameter definitions, see Maccio et al. (2007) 



Environment, -Di,o.i The environment in which the halo 
is residing is defined as the distance to the nearest FOF halo 
with mass greater than 0.1 M200, divided by the neighbour's 
Rzoo- This measure of environment is related to the inverse 
of the tidal force exerted by the neighbour and, contrary to 
other environmental parameters used in the literature, does 
not correlate with halo mass ( |Haas et aL|2011[ ). 



3 SIMULATIONS 

We use a sample of five cosmological N-body simulations 
with box sizes increasing by factors of two from L — 
25 h~ x Mpc to 400 ft -1 Mpc. The simulations were run with 
the TreePM code gadget-3 (|Springel|[2005[|. We adopted 



a set of cosmological parameters consistent with the 7-year 
WMAP data ( jKomatsu et al.||2009| ): (fi m , tt A , h, cr 8 , n s ) = 
(0.258,0.742,0.719,0.796,0.963). Glass-like cosmological 
initial conditions were generated at redshift z = 127 using 
the Zeldovich approximation and a transfer function gener- 
ated with CMBFAST (Seljak & Zaldarriaga 1996, v. 4.1). Each 
simulation contains 512 ,i collisionless particles. The particle 
mass is 1.16 x 10 7 M Q (L/25 h' 1 Mpc) 3 . The gravitational 
softening was set to l/25th of the mean inter-particle spac- 
ing, but was held fixed below z = 3. 

Catalogues of dark matter haloes were produced using a 
friends-of-friends (FOF) algorithm with dimensionless link- 
ing length b = 0.2. The FOF haloes were further analysed 
using the SUBFIND substructure finder (Springcl et al.|20uT |. 

Different parameters require different resolutions, with 
/sub and C200 being the most demanding. To ensure complete 
convergence for all properties, we only use haloes with at 
least 10 4 particles. We verified that using all haloes with at 
least 10 3 particles gives nearly identical results. The selected 
sample includes 1867 haloes (488, 534, 492, 306 and 47, re- 
spectively from the L = 25, 50, 100, 200 and 400 h' 1 Mpc 
simulations). 



4 RESULTS 

The relations between the different parameters and their 
relative importance are investigated using a non-parametric 
correlation analysis and a principal component analysis in 
§4.1 1 and §4.2[ respectively. 



4.1 Correlation analysis 

Figure [l] shows how the 9 halo properties correlate with each 
other. Each panel shows the distribution of haloes projected 
onto the 2 dimensions that are plotted along the axes. The 
colour scale indicates the logarithm of the fraction of haloes 
in a pixel. The value of the Spearman rank correlation coef- 
ficient, Rs, is printed in each panel. Values of Rs run from 
— 1 (perfect anti-correlation) to 1 (perfect correlation), while 
Rs — means no correlation. For reference, for our sample 
of 1867 haloes, values \Rs\ > 0.06 reject the hypothesis that 
no correlation is present at greater than 99% confidence. For 
the panels along the diagonal the two dimensions are iden- 
tical and we plot instead histograms of the values for the 
corresponding parameter with the right vertical axes show- 
ing the number of haloes per bin. Note that the figure is 
symmetric with respect to the diagonal. 

The strongest correlation, with Rs = 0.65, is between 
age (20.5) and concentration (C200), with older haloes being 
more concentrated. This correlation is expected, as the con- 
centration is thought to reflect the time since the assembly 



of the inner part of the halo (e.g. Zhao et al. 2003[ ). Substruc- 
ture (/sub), mass (M200), sphericity (s), and relaxedness (r) 
are all also strongly correlated (\Rs\ > 0.46) with concentra- 
tion. More concentrated haloes tend to have less substruc- 
ture, be less massive, more spherical, and more relaxed. This 
can also be understood, as more concentrated haloes are 
older and have thus had more time to relax and to tidally 
strip substructures. Moreover, tidal forces, and thus the rate 
at which substructures are destroyed, will be larger in more 
concentrated haloes. More concentrated haloes also tend to 
have lower spins (A) and to be more oblate (smaller T), but 
these correlations are less strong (_Rs ~ 0.3). Age shows 
nearly the same correlations as concentration, although it 
is interesting to note the anti-correlations with substructure 
and spin are much weaker for age than for concentration. 

The two strongest correlations that involve neither con- 
centration nor age are between substructure and spin (Rs = 
0.46) and between sphericity and mass (Rs = —0.40). The 
former correlation may reflect the fact that orbital motions 
of satellites contribute to the spin. Triaxiality and environ- 
ment appear to be the most independent parameters, corre- 
lating most strongly with sphericity (Rs — —0.33 for triax- 
iality and -0.19 for environment). 

Note that our measure of environment was designed to 
be insensitive to mass, which it indeed is (Rs = 0.07). It is 
thus quite different in nature from the environmental param- 
eters used in the literature, which all correlate strongly with 
mass ( Haas et al.|201l' l. The parameter -Di,o.i measures the 
distance to the nearest halo that is at least 0.1 times as mas- 
sive and expresses this distance in units of the neighbour's 
virial radius. Typical values are ~ 10 which, given that most 
neighbours have masses near the lowest allowed values, im- 
plies that the nearest neighbours are typically at distances 
of a few virial radii. The absence of strong correlations with 
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Figure 1. Correlations between different halo properties. The colour scale indicates the logarithm (base 10) of the fraction of haloes in 
the pixel. The Spearman rank correlation coefficient is shown in the top left of each panel. To guide the eye, values \R,s\ > 0.3 are printed 
red. The panels on the diagonal show histograms of the parameter values with the right vertical axes showing the number of haloes in 
the bins. Note that the figure is symmetric with respect to the diagonal. 



Difi.i may indicate that all aspects of the environment that 
do influence the properties of haloes are correlated with halo 
mass. 

It is interesting that all halo properties correlate more 
strongly with concentration than with mass (except for en- 
vironment, which, however, does not correlate much with 
either parameter). The difference is particularly striking for 
substructure and spin. While we find no evidence for a de- 
pendence of spin on mass (Rs = —0.04), spin is substan- 
tially anti-correlated with concentration (Rs = —0.32), as 



also found by Maccio et al. (20071. In agreement with Gao 



et al. (20041, substructure is nearly independent of mass 



(Rs = 0.09), but strongly anti-correlated with concentra- 
tion (Rs = —0.57). It thus appears that concentration, and 
the very closely related parameter age, may in some sense be 
more fundamental than mass. A better way to judge which 
(combinations of) parameters are most fundamental, is to 
carry out a PCA, which we will do next. 



4.2 Principal Component Analysis 

PCA is a statistical technique to find the number of inde- 
pendent parameters which are needed to account for the 
variance in a data set. Most of the variance in the data 
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Figure 2. Correlations between different halo properties and principal components. The colour scale indicates the logarithm (base 10) of 
the fraction of haloes in the pixel. The Spearman rank correlation coefficient is printed in the top left of each panel. For each parameter, 
the value of the strongest correlation is printed in red. The eigenvalue of each PC is shown along the horizontal axis. 



set can be represented by a subset of all eigenvectors of 
the covariance matrix. These eigenvectors are known as the 
principal components (i.e., PCI, PC2, etc.). The eigenvalues 
corresponding to the eigenvectors indicate whether or not 
that PC is important. If the variables are standardized (i.e. 
transformed to have mean zero and variance one), then the 
fraction of the total variance in the data set due to a given 
PC is equal to its eigenvalue divided by the total number of 
parameters. 

Except for / su b, which can be zero, we use the logarithm 
of the parameters. Before running the PCA, we standardize 
the halo properties by subtracting the means and dividing 
by the standard deviations. The results are insensitive to 



whether we take the logarithm or not, but standardizing is 
crucial because otherwise the result would depend on the 
choice of mass unit (note that all other parameters are di- 
mensionless) . 

Figure[2]shows the results of the PCA. Each panel shows 
the correlation between the halo property plotted along the 
vertical axis and the PC plotted along the horizontal axis. 
As in Figure [l] the values of the Spearman rank correla- 
tion coefficients are shown and the colour scale indicates the 
distribution of haloes. The eigenvalue of each PC is shown 
along the horizontal axis. 

PCI clearly stands out as most important, accounting 
for no less than 37% of the total variance, compared to 15% 
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for PC2. PCs 3 and 4 also account for > 10% of the variance, 
but the remaining PCs are not that significant. However, to 
account for 90% of the total variance, we require no less than 
6 PCs, a suprisingly large number. Even the least important 
PC still accounts for more than 2% of the variance. Clearly, 
dark matter haloes are complicated objects whose structure 
cannot be described using a small number of parameters. 

PCI correlates very strongly with concentration (Rs = 
—0.86), which indicates that this parameter is most funda- 
mental. Environment, spin, and particularly triaxiality also 
correlate significantly with PCI, but they correlate more 
strongly with other PCs. All other parameters correlate most 
strongly with PCI. This is consistent with the fact that its 
eigenvalue is by far the greatest. This confirms what Figure[l] 
suggested: of the parameters that we consider, concentration 
is most fundamental. 

PC2 correlates very strongly with spin (Rs = 0.72), 
strongly with mass (Rs = —0.56), but not at all with con- 
centration. PC3 only correlates strongly with environment 
(Rs = 0.68). PC4 correlates most strongly with triaxiality 
(Rs = 0.64), but is also strongly correlated with environ- 
ment (R s = -0.55). 

Hence, concentration, age, relaxedness, mass, and sub- 
structure all belong to a single family of parameters, which 
however, still contains a large amount of scatter. PCs 2-4 
together account for the variance that is not linked with 
concentration. That part of the variance is mostly due to 
the scatter in spin, environment, and triaxiality, the three 
parameters that are most independent of the main family of 
halo properties. 

As this paper was in the final stages of preparation, 



Skibba & Maccio (20111 posted a preprint of an indepen- 
dent, related study. While they also carried out a PCA, 
albeit with a somewhat different set of parameters, they 
did not present a correlation study. In agreement with our 
PCA results, they find that concentration is more funda- 
mental than mass. However, they also find that relaxedness 
is about as important as concentration, whereas we find that 
the variation with relaxedness that does not trace concen- 
tration only features strongly in PC5, which accounts for 
only 5% of the total variance. This difference may be due 
to the fact that their PCA includes multiple measures of 
relaxedness and that they did not consider age, which is 
closely related to concentration. Another cause of discrep- 
ancy is the use of different definitions for environment. They 
use the overdensity in a fixed aperture (of 8 Mpc/h), which 
has been shown to correlate very strongly with mass ( Haas 
et al.|20TT l, whereas our environmental parameter is insen- 
sitive to mass. 

To conclude, both the correlation analysis and the PCA 
demonstrate that concentration, age, substructure, mass, 
sphericity, and relaxedness are closely related, with con- 
centration being most fundamental. Triaxiality, spin, and 
(mass-independent) environment are more independent, al- 
though spin correlates strongly with substructure and both 
spin and triaxiality are substantially correlated with concen- 
tration. While the scale of a halo is set by its mass, all other 
properties are more closely related to concentration. 
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